Two-Particle Self-Consistent Approach to Anisotropic Superconductivity 
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A nonperturbative approach to anisotropic superconductivity is developed based on the idea of 
two-particle self-consistent (TPSC) theory by Vilk and Tremblay. A sum-rule which the momentum- 
dependent pairing susceptibility satisfies is derived. An effective pairing interaction between quasi- 
particles is determined so that the susceptibility should fulfill this exact sum-rule, in which fluctu- 
ations belonging to different symmetries couple at finite momentum. It is demonstrated that the 
mode coupling between d-wave and s-wave pairing fluctuations leads to suppression of the d-wave 
fluctuation near the Mott insulator. 

PACS numbers: 74.20.-z, 71.10.Fd 



I. INTRODUCTION 

The anisotropic superconductivity originating in the 
Coulomb repulsion has been a topic of continuous inter- 
est in strongly correlated electron systems. In the typical 
situation described by the Hubbard model, the super- 
conducting phase lies in the regime with two comparable 
energy scales: the bandwidth W and the Coulomb re- 
pulsion £/[l[. This intermediate regime is unreachable 
by the Monte Carlo simulations on account of the severe 
sign problem Q. Therefore, addressing the anisotropic 
superconductivity has been a challenging issue both in 
numerically and analytically. 

The elemental approach is the perturbation theory 
with respect to U . The random-phase approximation 
(RPA) gives an intuitive picture of the effective interac- 
tion arising from A more systematic approxi- 
mation is the fluctuation exchange (FLEX) approxima- 
tion, in which all physical quantities are derived from 
the Luttinger-Ward functional 0, These theories give 
us basic information on the dominant pairing fluctuation 
in models and realistic materials [1]. In addressing the 
strong-correlation regime, U > W , however, the pertur- 
bative treatment of U is not sufficient. As a result, the 
superconducting fluctuation survives close to the Mott 
insulator [l], The local correlation, which is incorpo- 
rated in the self energy, is responsible for the suppression 
of the superconducting fluctuation. 

On the other hand, the strong local correlation giving 
rise to the Mott transition can be treated by the dy- 
namical mean- field theory (DMFT) 0, [ljj . Its cluster ex- 
tensions enable us to address the momentum-dependent 
vertex part which is indispensable for the anisotropic 
superconductivity (Til [T3 |. To exclude the size effect, 
it is necessary to take larger clusters than the minimal 
sizepjj]. Other kinds of extensions without clusters have 
also been proposed to incorporate the influence of the 
long-range spatial fluctuations (T^ - [2fjj . 

A non-perturbative approach referred to as two- 
particle self-consistent theory (TPSC) has been devel- 
oped by Vilk and Tremblay (U [Hj]. This theory starts 



from the two-particle fluctuation in contrast to the above 
approaches, where the self-energy is evaluated prior to 
the two-particle fluctuation. In this sense, TPSC inher- 
its the idea of the phenomenological method named self- 
consistent renormalization (SCR) theory, which has suc- 
ceeded in describing the quantum critical phenomena [23l— 
[25j . TPSC consists of two parts. First, an effective quasi- 
particle interaction is derived from the double occupancy 
so that the spin and charge susceptibilities satisfy the 
exact sum-rule. In this sum-rule, the local correlation is 
taken into account by the quantity {n,jin a i), which corre- 
sponds to the double occupancy for different spins and is 
reduced to (n a } for the same spins by the Pauli princi- 
ple. The interaction is assumed to be independent of the 
energy and momentum but dependent on the spin com- 
ponents. With this spin-dependence and the sum-rule, 
they have achieved a reasonable description of the spin 
and charge fluctuations taking account of the local corre- 
lation. In the next step, the two-particle susceptibilities 
thus obtained are used to derive the self-energy, which 
therefore incorporates influence of the collective modes. 

After the development, TPSC has been applied to su- 
perconductivity. The attractive Hubbard model has been 
investigated to discuss the isotropic pairing [26l |27|. The 
superconductivity arising from the repulsive interaction 
has been discussed by Kyung et al. [28| The superconduct- 
ing transition temperature and the single-particle prop- 
erties have been discussed. In this paper, we address the 
anisotropic superconductivity in a way different from ref. 
(28j . Namely, we extend the first step of TPSC (effective 
interaction) to the superconducting fluctuations. To this 
end, an exact sum-rule which the pairing susceptibility 
satisfies is derived. Using this sum-rule, the pairing inter- 
action as well as the pairing susceptibility is determined. 

To make the discussion concrete, we consider in this 
paper the two-dimensional Hubbard model 
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where e/~ = — 2t (cos k x + cos k y ) — At' cos k x cos k y , and we 
take t = 1 as a unit of energy. The system size is N = L 2 



2 



and we adopt the periodic boundary condition. 

This paper is organized as follows. We begin with a 
review of TPSC in the next section. Section ITTT1 describes 
its extension to superconductivity. We first derive a sum- 
rule and then an equation for an effective pairing inter- 
action is proposed. Numerical results are given in Sec- 
tion IIV1 We summarize in Section [V] with some discus- 
sions. 



This relation expresses the fact that the effective inter- 
action is reduced in accordance with the decrease of the 
probability that two electrons having opposite spins exist 
at the same site. Thus equations are closed, and t/ c h, U sp 
and (n-|-nj.) are evaluated self-consistently for given n and 
U. We note that Eq. (|9]) breaks the particle-hole symme- 
try and should be used only for n < 1. The range n > 1 is 
considered through the particle-hole transformation into 
n < 1. 



II. TWO-PARTICLE SELF-CONSISTENT 
(TPSC) THEORY 



III. EXTENSION TO SUPERCONDUCTIVITY 



In this section, we review TPSC by Vilk and 
TremblaydlHi]. 

We define the susceptibility by 



(2) 



The spin and charge susceptibilities are defined from Xaa > 
respectively, by 

Xs P = 2(xtf -Xf4-)> Xch = 2(xtt + Xn)- ( 3 ) 

The TPSC takes account of the spin-dependence of the 
effective interaction. We express this coupling constant 
by U a<T '. We assume = U~n and U\\, = and 
introduce 



U, 



sp 



U n - U- 



u ch = u n + u. 



(4) 



With this effective interaction, the susceptibility in the 
RPA is given by 



x^ psc) (q) = 



2X0 (q) 



1 - U sp Xo(q)' 
2X0 (g) 

i + u ch xo(q) ' 



(5) 
(6) 



where X o(q) = -(T/N)J2 k G(k)G(k + q) and G(fc) = 
l/(ie n + \i — £fc). We have introduced the notations 
k = (k,ie n ) and q — (q,iv n ) with e n and v n being 
the fermionic and bosonic Matsubara frequencies, respec- 
tively. The effective coupling constants are determined so 

that Xsp PSC ^(<z) and Xa? (q) satisfy the exact sum-rule. 
The self-consistency equations thus read 



T 



^E4 PSC) (9)-™-2(n t rH), 



x ch PSC) (q) = n + 2 ( n t n i) - n 2 , 



(7) 
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where n is the particle number per site. We have assumed 
( Tl t) = ( n i) = n /2 an d have used the relation (n^) = 
(n£), which expresses the Pauli principle. Provided that 
the double occupancy (n^n±) is given, U ap and C/ C h are 
determined from the above equations. In the TPSC, the 
following assumption that connects (n-j-n^) with U sp is 
invoked 



(«t)W) 



u. 



(9) 



A. Sum-rule for Pairing Susceptibility 

In the TPSC reviewed in the previous section, the sum- 
rule for the susceptibility plays the main role. To apply 
this idea to the anisotropic superconductivity, we first 
derive a corresponding sum-rule for the superconducting 
fluctuations. 

We consider the Cooper pair with the total momen- 
tum q. Although the present interest is on the zero- 
momentum pairing, finite momentum is necessary to con- 
struct a sum-rule which the pairing susceptibility satis- 
fies. The annihilation operator B aq for the symmetry a 
is defined by 

fe 

= W^e-^-(^+«A (10) 

ij 

Here g a (k) denotes the orbital part of the pair electrons, 
and is classified in terms of the irreducible representa- 
tions of the point group. (g a )ji is the Fourier trans- 
form of g a {k): {g a ) Jt = N- 1 ^ k g a (k)e ih < a i- B *\ We 
note that the range of the vector q is q x = (ATr/L)n x 
with n x — 0, ••• ,L — 1, since the periodicity of the 
center-of-mass coordinate (Uj + Rj)/2 is half of iJj. 
Hence, for anisotropic pairing in general, B q+ a ^ B q 
and B q+2G = B q . 

The pairing dynamical susceptibility is defined by 

Poof {q) = f dr^(B aq {r)Bl q )e^ 



^Y.9«(k)P{k,k' ]q )g* a ,{k'), (11) 
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where 



P(k,k';q) 

dT(c k+ ^(T)c_ k+ii (T)cl k , + u cl +it )e l ^ T . 



(12) 



The off-diagonal components of P aa i are finite for q ^ 0. 
For q = 0, on the other hand, P aa ' is block-diagonal with 



TABLE I: The even-parity orbital function g a {k) and the 
corresponding real-space vector r a = Rj — R4 for which 

(fla)i« / 0. 



symmetry a 


g a (k) 


r a 


Aig 


1 


(0,0) 


A 2g 


\/2(sin 2fc x sin 2fej, — sin k x sin ) 


(2,1) 


Bi g 


cos k x — cos fcy 


(1,0) 


B 2g 


2 sin fc^ sin fc K 


(1,1) 
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DO G— O ,0, 



respect to the irreducible representations in the point 
group. 

The Fourier transform is defined by P(r,r) = 
(T '/N)Y Jq . n P{q)e lq ' r ~' WnT ', where the prime indicates 
that each element of q takes the values expressed by 
q x = (4ir/L)n x as mentioned above. Then the equal- 
time local component, P aa '(f = 0,r = —0), is expressed 
as 

ijlm 

= Qaa' j (1^) 

where 8 is Kronecker's delta extended to satisfy the pe- 
riodic boundary condition. Q aa ' is block-diagonal in a 
manner similar to P aa > (q — 0, iv n ). Eq. ([T3|) is the sum- 
rule which the momentum-dependent pairing susceptibil- 
ity satisfies. 

In the square lattice, g a is classified by the irreducible 
representation of the point group D4h- Table U shows 
the even parity representations with the smallest r a in 
each irreducible representation. Here r a denotes one of 
vectors r a = Rj — Ri for which {g a )ji ^ 0. We have 
chosen the normalized condition TV -1 ^ fe |g Q (fc)| 2 = 1. 
Explicit expressions for Eq. (|13p are then given by 

QAig.Aig = (n>t n l)> ( 14 ) 

QBlg,Blg = («t(*) n i) - ( cr +(*) (7 -), ( 15 ) 

QB2g,B2 g = (rif(x + - (<r+(x + y)a-) 

-2(c\(x + y)c r (x)clc±{y)}, (16) 

where x and y denote the primitive translation vectors. 
Corresponding diagrams are shown in Fig. [T] For r a ^ 
0, Q aa consists of two-site terms and four-site terms. 
When r a — (2n + 1,0), the four-site term is absent as in 
Eq. (|15|) . The two-site terms may be rewritten in terms 
of the spin and charge susceptibilities as 

(n t (r„)n;) - (<7 + (r Q )cr_) 

= \ [Xch(r a , t = 0) - 3 X s P (r Q , r = 0) + n 2 } . (17) 
Here we have assumed isotropy in the spin space. 



FIG. 1: Diagrammatic representations of Q aa for a = Ai g , 
Big and B2 g . 

The explicit forms of Qaa 1 in Eqs. (|14l) - (fTH)) confirm 
the importance of the sum-rule (fT3|) . The double occu- 
pancy (n^ni), appeared in Eq. (fll)) , is one of the most 
important quantities in systems with the local repulsion. 
For anisotropic pairings, on the other hand, Q a a con- 
sists of spin and charge correlations as the effective pair- 
ing interaction in the RPA0. Therefore, Q aa reflects 
the influence of either the local correlation or the spin 
correlation depending on the basis a. 

B. Self-Consistent Equations 

We express the pairing susceptibility P aa ' (?) in a RPA 
form phenomenologically. To this end, we consider the 
following effective interaction 

Hint = £ V ( k ~ fc ') C l+a t c -fc+a; c -fc' + f 4. c fc' + § f 
kh' q 

(18) 

Here we assume that the coupling constant depends only 
on the momentum transfer k — k' . Then, V can be writ- 
ten in terms of the irreducible representation g a {k) as 
follows: 

7(fc-fc0 =£«(%<*(*')■ (19) 

a 

Here we have assumed a single basis for each irreducible 
representation. Actually, V may have off-diagonal el- 
ements between (/ Q 's belonging to the same irreducible 
representation. 

In the RPA, the two-particle Green function P(k 1 k': q) 
in Eq. (fT2j) is given by 

P^ s \k,k';q) = P (k;q)6 kk , 

~ Jj E P °( fc ; ^( fc - k "W, k'- q), (20) 
fc" 

where 

P (k; q) =T^G(fc+ |,ie ro )G(-fe + |, -ie m + iv n ). 

m 

(21) 
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Inserting Eqs. (fT9j and (J20]) into Eq. ([TT]), we obtain the 

expression for P^,\q) in the RPA. To make the notation 
simple, we use a matrix form with respect to a indices 
and denote the matrix by a hat. Then P^ eS \q) is ex- 
pressed as 



P {eS) {q) = [Po{q)- 1 + V] 



(22) 



where (V) aa > — S aa 'V a and Po, a a'{q) is defined from 
Po(k;q) in a manner similar to Eq. (fTT|) . We substitute 
p( eff )(g) i n to the exact sum- rule, Eq. (TlU)) . The equation 
reads 



T 
N 



Y^P {cS \q)e w " 0+ =Q. 



(23) 



This equation determines the effective pairing interac- 
tions V a from Q aa ' , which consists of double occupancy 
(n^ni) or equal-time correlations such as Xsp( r a,i~ = 0) 
depending on the symmetry a. 

To see the tendency of the solution, let us consider 
the simplest situation where the off-diagonal components 
of Po(q) are neglected. Under this approximation, if 
Q aa is larger than the non-interacting value, the effec- 
tive interaction V a may be attractive to enhance Paa (q) 
in Eq. (|23|) . For Ai g symmetry, for example, the re- 
pulsive U reduces <5Ai g ,Ai g in Eq. (TT4"]) . leading to a 
repulsive effective interaction. For Bi g symmetry, on 
the other hand, the antiferromagnetic spin fluctuation of 
q = (7r, 7r) around the half-filling gives a negative value 
for Xs P (£, t = 0) to increase Qriig.Big in Eq. (fT5|) . Hence 
the pairs with the Bi g symmetry is enhanced as expected. 

The issue of interest is thus how the solution is affected 
by the off-diagonal components of Po(q), which has fi- 
nite values even between different symmetries at q ^ 0. 
The off-diagonal susceptibility couples the pairing fluctu- 
ations in different symmetries and therefore, the Ai g and 
Bi g fluctuations are not independent any more. Numer- 
ical solutions are given in the next section with central 
attention on this point. 



using Eq. (fTTj) . we obtain the following expression for 
Eqs. (HU) and ([T^]): 



)(2-site) 



W (r tt ,r = 0)-3 X W(r a ,r = 0) 



Xch 



+ xo(r a ,r = 0) + L^' P .a a (q)e^ 0+ . (24) 



The procedure for solving Eq. (123[) is summarized as 
follows. We first solve TPSC equations, Eqs. (El)-®, 
to obtain (n-|-nj,), U sp and U c h- In the next step, these 
results and Eq. are used as the input to Eq. (|2"3"j) . 
thus obtaining V and P^ eS ^(q) for arbitrary q — (q,iv n ). 



D. Mermin- Wagner Theorem 

We conclude this section by demonstrating that 
the present approach satisfies the Mermin- Wagner 
theorem [H H| similarly to the TPSC. The proof in 
ref. (2lj . which concerns the magnetism, is also appli- 
cable to superconductivity. To see this, we observe that 
Qaa in the right-hand side of Eq. (|23|) is a quantity of 
the order of unity. On the other hand, the left-hand side 
of Eq. (f2"3"]l diverges toward the second-order phase tran- 
sition in one- and two-dimensions. This follows from the 
assumption that the static component of p( cff ) (q) takes 
the form 



p(eff) (g . 



<7o 0) 



1 



(25) 



near the critical point, where £ denotes the correlation 
length. Therefore the superconducting phase transi- 
tion of the second order is forbidden in one- and two- 
dimensions provided the sum-rule (fT3f is satisfied. This 
property will be confirmed numerically in the next sec- 
tion. 



IV. NUMERICAL RESULTS 



C. Approximation for Equal-Time Correlations 



A. Equal-Time Correlations 



In the present framework, the key quantity is Q a a', 
which consists of equal-time short-range correlations. 
Provided that Q a a' is given, the effective pairing inter- 
actions V a are determined from Eq. (|23)) . To evaluate 
Qaa 1 , we may use external numerical methods such as 
the exact diagonalization and the quantum Monte Carlo. 
Instead, we here take advantage of TPSC results reviewed 
in Section [TT] so as to make the equation closed. 

As shown in Eqs. (TTS")) and p^|) . Q aa for the anisotropic 
pairing consists of two-site correlations and four-site cor- 
relations. For the two-site correlations, we use the 
TPSC results, Xsp psc ' ) and xih PS °^ m Section |TTJ Concern- 
ing the four-site term, we replace them by their non- 
interacting values. Eliminating the four-site terms and 



In this section, we show numerical results for supercon- 
ducting susceptibilities. We use, as the input to Eq. (|2"3"1) . 
TPSC results in Section |TTJ First, we check the accuracy 
of this input by comparing with the exact diagonalization 
(ED) method. 

Fig. [5] shows Qaa for a — Ai g and Bi g as a function 
of the particle number n per site. The system size L is 
L = 128 for TPSC and L = 4 for ED. The tempera- 
ture dependences of TPSC results differ in the symme- 
tries: correlations for Ai g symmetry is suppressed with 
decreasing temperature, while that for Bi g symmetry is 
enhanced. These tendencies are consistent with the sit- 
uation caused by the repulsive U as discussed in Sec- 
tion IIII Bl Among three values of T shown in Fig. [21 
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FIG. 2: (Color online) Q aa evaluated in TPSC (Eq. (23)1) as 
a function of n for t = 0, U = 8 and L = 128. The dots are 
results computed by ED with L = 4. 



FIG. 3: (Color online) The effective pairing interaction V a at 
T = 0.1 for t' = 0, U = 8 and L = 128. The arrow indicates 
a minimum of Vkig in approximation (ii) . 



T = 0.1 is the closest to the ED results. This corre- 
spondence is reasonable in view of the fact that, in the 
system with L = 4, the energy levels of the free electrons 
have gaps of the magnitude 0.25. 

B. Static Susceptibilities 

We proceed to the results for the effective pairing in- 
teractions and the pairing susceptibilities. The solution 
of Eq. (|23l) depends on the choice of g a . We apply the 
following two approximations: 

(i) A lg or Bi g : We neglect the off-diagonal compo- 
nent of Pq. In this case, the matrix equation (|23|) 
becomes diagonal. 

(ii) A\g and B\ g : We solve the 2 by 2 matrix equation 
with a — Aig and Bi g . 

The approximation (ii) is the minimal choice to in- 
clude both the local correlation and the spin fluctuation. 
Comparison of results in the above two approximations 
makes it clear the influence of the off-diagonal component 
Po,aa' i which couples pairing fluctuations in different 
symmetries, in the present equation. We shall not present 
results including other symmetries and long-range pair- 
ings, since they make only a quantitative change. 

Fig. [3] shows the effective pairing interactions V a as a 
function of n. The result in the approximation (i) agrees 
with the analysis in Section MI Bl Vai s is repulsive be- 
cause of suppression of (n-^nj,), while Vei g is attractive 
due to the enhancement of the antifcrromagnetic correla- 
tion between spins on the nearest neighbor sites. The at- 
traction in this approximation, however, turns out to be 
inaccurately much enhanced so that the superconducting 
fluctuation predominates the antiferromagnetic fluctua- 
tion (this situation becomes clearer in Fig. [3]). Inclusion 




0.5 0.6 0.7 0.8 0.9 1 



n 

FIG. 4: (Color online) The inverse of the static susceptibility 
Paa (0,0) with a = Big as a function n for t' = 0, U = 8 
and L = 128. 



of the off-diagonal component of Pq suppresses the at- 
traction in Big symmetry. We can see, in approximation 
(ii), that Veig turns to repulsive, the extent of which 
is conspicuous around n = 1. Nevertheless, T^Big still 
has tendency to attraction, making a minimum around 
n = 0.9 as indicated by the arrow in Fig. [3] 

With the effective interaction V a shown above, 

Paa'\ ( l) f° r arbitrary q and v n can be evaluated. We 
here show results only for q — and v n — in Bi g 
symmetry and discuss the tendency to the superconduc- 
tivity. The inverse of Pi C a\o, 0) with a = Big is shown 
in Fig. |4j In the approximation (i) , the susceptibility is 
nearly diverged in a wide range of n > 0.7 at T = 0.1. 
In actual, it does not diverges exactly, since the present 
theory satisfies the Mermin- Wagner theorem. Compared 
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L= 128 



L=128 
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0.7 



0.8 



0.9 
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FIG. 5: (Color online) A phase diagram of the "quasi-two 
dimensional" system calculated in the present theory (Bi g 
superconductivity) and in TPSC (magnetism) for U — 8. AF- 
C indicates the antiferromagnetism with q = (ir, ir) and AF- 
IC indicates ordering with the vector other than (ir,ir) The 
horizontal dotted line indicates the lowest temperature T = 
0.01 in the calculation. The system size is L — 128 (closed 
symbol) and L — 64 (open symbol), (a) t' = and (b) 
t' = -0.1. 



to this result, the susceptibility in the approximation (ii) 
is strongly suppressed, and its inverse first touches the 
zero at T ~ 0.01. It is remarkable that, in contrast to 
the approximation (i), the Bi g fluctuation does not go to- 
ward divergence around n = 1, accompanying the strong 
suppression of the A lg fluctuation. 

Although there is no phase transition in the present 
calculation in two dimensions, to see the tendency of the 
fluctuations, we define a "transition temperature" as the 
temperature at which the susceptibility equals 10 3 . Wc 
refer to the resultant phase diagram as "quasi-two di- 
mensional" phase diagram hereafter. The phase diagram 
for t' = is shown in Fig. [5ja). The antiferromagnetic 
transition temperature Taf calculated in TPSC is also 
plotted. The superconducting transition temperature T c 
in the approximation (i) is comparable to Taf- In the ap- 



proximation (ii), on the other hand, the superconducting 
phase appears only around n — 0.9 as is expected from 
Veig in Fig. 2J However, it turns out that T c is always 
smaller than Taf- This apparently unreasonable result is 
ascribed to the separate evaluation of Taf and T c , which 
will be discussed in the next section. We also show the 
"quasi-two dimensional phase diagram" for t' = —0.1, 
which is expected to be a typical situation in the cuprate 
superconductors [l|. The range n > 1 has been evaluated 
by transforming into n < 1. The cusp at n = 1 is due to 
the use of Eq. ©, which breaks the particle-hole sym- 
metry. We can see an enlargement of the superconduct- 
ing phase in hole-doped regime while it is comparable to 
or even narrower than that of t' = in electron-doped 
regime. However, the problem that T c does not overcome 
Taf still remains in this parameter. 



V. SUMMARY AND DISCUSSIONS 

We have extended TPSC to anisotropic superconduc- 
tivity. The sum-rule, Eq. (1 13[) , is the key equation, which 
relates the superconducting susceptibility to the spin and 
charge correlations expressed by Q aa > . The effective pair- 
ing interactions are determined so that the RPA-type 
phenomenological susceptibility satisfies the exact sum- 
rule. We remark that, in the present matrix equation, 
the off-diagonal susceptibility Po, aa '{q) plays an impor- 
tant role in coupling the superconducting fluctuations in 
different symmetries at finite q. As a result of this "mode 
coupling" , the d-wave pairing fluctuation is suppressed 
around half-filling accompanying the s-wave pairing fluc- 
tuation, which is strongly suppressed by the local repul- 
sion U. 

Describing the suppression of the d-wave fluctuation 
near half-filling has been one of issues in the perturba- 
tion theories. An essential factor lacking there is the 
local correlation effect, which the self-energy correction 
is one of factors responsible for. In the present equa- 
tions, on the other hand, the suppression is achieved by 
the coupling with the s-wave fluctuation, which directly 
reflects the local repulsion U. In other word, the lo- 
cal correlation effect is incorporated within two-particle 
quantities through the mode coupling between supercon- 
ducting fluctuations in different symmetries. 

Although the <i-wave fluctuation is indeed suppressed 
around half-filling, T c seems inaccurately small compared 
to Taf- This difference is ascribed to the structure of 
the present equations. In the present framework, Xsp 
and x c h are first calculated to obtain Taf and with use 
of them, the pairing fluctuations are evaluated. Hence, 
the feedback from the superconducting fluctuation to the 
spin fluctuation is not taken into account. This consid- 
eration brings us a possible improvement of the present 
theory: to treat the spin and superconducting fluctua- 
tion equally, an anisotropic vertex part giving rise to the 
unconventional spin-density wave [31] should be included, 
and we may further bring the parquet formalism to treat 
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both the fluctuations on an equal footing [32l - [34j 

Finally, we again remark the importance of the sum- 
rule Eq. fp~3|) for pairing susceptibility. Theories which 
satisfies this sum-rule follows the Mermin- Wagner theo- 
rem. The perturbation theories such as RPA and FLEX 
as well as DMFT and its cluster-extensions do not satisfy 
this sum-rule. This sum-rule could be one of directions 
in developing a framework addressing the anisotropic su- 
perconductivities. 



Acknowledgments 

The author would like to thank H. Kusunose for useful 
discussions and comments on the manuscript. He is also 
indebted to J. Schmalian, T. Takimoto, D. Vollhardt and 
H. Yokoyama for stimulating discussions and J. Nasu for 
a technical advice. The author is supported by JSPS 
Postdoctoral Fellowships for Research Abroad. 



For a review, see Y. Yanase, T. Jujo, T. Nomura, H. [18 
Ikeda, T. Hotta, and K. Yamada, Phys. Rep. 387, 1 
(2003). [19 
E. Y. Loll Jr., J. E. Gubernatis, R. T. Scalettar, S. R. 
White, D. J. Scalapino, and R. L. Sugar, Phys. Rev. B [20 
41, 9301 (1990). 

N. F. Berk and J. R. Schrieffer, Phys. Rev. Lett. 17, 433 
(1966). [21 
P. W. Anderson and W. F. Brinkman, Phys. Rev. Lett. 
30, 1108 (1973). [22 
D.J. Scalapino, E. Loh, Jr., and J.E. Hirsch, Phys. Rev. 
B 34, 8190 (1986). 
N. E. Bickers and D. J. Scalapino, Ann. Phys. (NY.) [23 
193, 206 (1989). 

N. E. Bickers and S. R. White, Phys. Rev. B 43, 8044 [24 
(1991). [25 
T. Takimoto and T. Moriya, Phys. Rev. B 66, 134516 
(2002). [26 
W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989). [27 
A. Georges, G. Kotliar, W. Krauth and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). [28 
T. Maier, M. Jarrell, T. Pruschke, and M.H. Hettler, Rev. 
Mod. Phys. 77, 1027 (2005). [29 
M. Potthoff, Eur. Phys. J. B 32 (2003) 429; M. Pot- 
thoff, M. Aichhorn, and C. Dahnken, Phys. Rev. Lett. [30 
91 (2003) 206402. [31 
T.A. Maier, M. Jarrell, T.C. Schulthess, P.R.C. Kent, 
and J.B. White, Phys. Rev. Lett. 95, 237001 (2005) [32 
H. Kusunose, J. Phys. Soc. Jpn. 75 (2006) 054713. 
A. Toschi, A. A. Katanin, and K. Held, Phys. Rev. B 75, 
045118 (2007). [33 
K. Held, A. A. Katanin and A. Toschi, Prog. Theor. Phys. 
Suppl. 176, 117 (2008). 
A. A. Katanin, A. Toschi, and K. Held, Phys. Rev. B 80, [34 
075104 (2009). 



C. Slezak, M. Jarrell, Th. Maier, and J. Deisz, J. Phys.: 
Condens. Matter 21 435604 (2009). 

A. N. Rubtsov, M. I. Katsnelson, A. I. Lichtenstein, and 
A. Georges, Phys. Rev. B 79, 045133 (2009). 
H. Hafermann, G. Li, A. N. Rubtsov, M. I. Katsnelson, 
A. I. Lichtenstein, and H. Monien, Phys. Rev. Lett. 102, 
206401 (2009). 

Y. M. Vilk and A.-M. S. Tremblay, J. Phys. I (Paris) 7, 
1309 (1997). 

A. -M. S. Tremblay, in " Theoretical methods for Strongly 
Correlated Systems", ed. by A. Avella and F. Mancini, 
(Springer Verlag, 2011) arXiv:1107.1534| . 

T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 
netism (Springer, Heidelberg, 1985). 
T. Moriya and K. Ueda, Adv. Phys. 49, 555 (2000). 
H. Kondo and T. Moriya, J. Phys. Soc. Jpn. 78, 013704 
(2009). 

S. Allen and A.-M. S. Tremblay, Phys. Rev. B 64, 075115 
(2001). 

B. Kyung, S. Allen, and A.-M. S. Tremblay, Phys. Rev. 
B 64, 075116 (2001). 

B. Kyung, J.-S. Landry, and A.-M. S. Tremblay, Phys. 
Rev. B 68, 174502 (2003). 

N.D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 

S. Coleman, Commun. math. Phys. 31, 259 (1973). 

H. Ikeda and Y. Ohashi, Phys. Rev. Lett. 81, 3723 

(1998). 

C. de Dominicis and P. C. Martin, J. Math. Phys. 5, 14 
(1964); C. de Dominicis and P. C. Martin, J. Math. Phys. 
5, 31 (1964). 

B. Roulet, J. Gavoret, and P. Nozieres, Phys. Rev. 178, 
1072 (1969); P. Nozieres, J. Gavoret, and B. Roulet, 
Phys. Rev. 178, 1084 (1969). 

H. Kusunose, J. Phys. Soc. Jpn. 79, 094707 (2010). 



